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Abstract 

The quantum properties of solitons at one loop can be related to phase shifts 
of waves on the soliton background. These can be combined with heat kernel 
methods to calculate various parameters. The vacuum energy of a CP(1) 
soliton in 2 + 1 dimensions is calculated as an example. 
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I. INTRODUCTION 



Solitons arise in many field theories and their particle-like behaviour has prompted dis- 
cussions of their quantum properties 0. Recently, some old results of Schwinger 0] have 
been developed into a useful numerical scheme for calculating the vacuum energy of a soliton 
at one loop order H^]. The aim of this paper is to develop these ideas further using heat 
kernel methods and discuss some applications. 

The one loop corrections to the soliton energies are likely to be most significant in 
situations where different classical soliton solutions have the same energy This happens 
in an important class of field theories, where the classical solutions fall into topologically 
seperate families and the solutions in each family saturates an energy bound. The prototype 
for this behaviour was the BPS monopole solution B, but similar solutions play an important 
role in superstring theories || . 

The CP(1) model in 2 + 1 dimensions is one of the simplest models of this type ]7|||. 
The calculation of the vacuum energy of the single soliton solution provides an instructive 
example of the general technique. A useful comparison can also be made with the CP(1) 
model in two dimensions, where the one loop contribution to the path integral from instanton 
solutions has been obtained analytically [pUlQj. 



II. HEAT KERNEL METHOD 

The one loop quantum properties of a soliton solution depend on the normal mode 
frequencies of classical perturbations about the solution. It will prove convenient to impose 
boundary conditions at a fixed radius R on the perturbations in order to obtain a discrete 
spectrum and then take the limit R —> oo. 

The heat kernel is defined by 

K(t)=Y,e- k2 \ (1) 

where the sum extends over the discrete spectrum with values k 2 . For large R, the normal 
modes with real k approach trigonometric functions of kR + 0, where is a constant phase 
depending on k and a set of other parameters 7. If no solitons are present, then the normal 
mode frequencies (denoted by a superscript zero) are given by 

k {0) R + <j) (0) =nn + (3, (2) 

where (3 is a phase that depends on the boundary conditions. In the presence of the soliton, 
the real frequencies are given by 

kR + = nvr + (3. (3) 

In the limit R — > 00, k — > k^°\ and the heat kernel can be expressed in terms of the phase 
shift by 

K{t) - {t) = ^-J^dk e~ k2t kt <y 7 (Jfe) (4) 



2 



where <5 7 = <fi - 

In fact, the n — 1 states can disappear from the continuum in the R — > oo limit if 
<5 7 (0) 7^ and the derivative 8' (0) = 0. In this case the phase shift must be displaced onto 
a new branch. Furthermore, bound states have their own asymptotic behaviour and need 
to be added as a seperate contribution to equation (||). If there are n bound states with 
k = 0, equation (|) becomes 

Q POO 

Kit) - K (0) {t) = n + - / dke- k2t kty]{5Jk) - 6 JO)). (5) 

7T JO 



Levinson's theorem |TTJ implies that the two expressions and (|^) are usually identical. 



However, in some situations (including the CP(1) soliton), Levinson's theorem does not 
apply and equation (|5|) in the one that must be used. 

An important feature of the heat kernel is the well-understood behaviour of the t — > 
limit in d dimensions ||T^|.[T^ . 



m-r^Y^Bne (6) 

n=0 

where the heat kernel coefficients B„ determine the one loop ultra-violet divergencies of the 
theory. Finite parts of physical quantities can be obtained from a regulated heat kernel, 

(d+l)/2 

K rcg (t) = K(t)-r d / 2 £ B n t n (7) 

n=0 

The heat kernel coefficients are integrals of local polynomials of the background fields and 



their derivatives. Explicit results are known for the first few coefficients p^JT5[| . 

The first term in the asymptotic expansion is equivalent to the free heat kernel K^°\ A 
simple comparison shows that inserting the first Born approximation to the phase shifts Sjp 
into equation (f|) produces next term in the asymptotic expansion with the coefficient B\. 
However, the second order Born approximation 8^ involves a double integral which contains 
contributions to B2 and higher terms. The part responsible for B2 can be isolated by taking 
a derivative expansion and keeping only the leading term. This allows us to express the 
regulated heat kernel in terms of phase shifts, 

K Teg (t) = n + -j^ dk kte~ k 1 X> 7 (A;) - £ 7 (0)) (8) 

where 5 7 = 8 1 —8^' subtracting the n'th order Born approximations to the phase shift 
keeping only local terms with m derivatives and In + 4m < d + 1. 

The heat kernel can be used to find the vacuum energy of a soliton in d + 1 dimensions 
or the effective action of an instanton in d dimensions. The vacuum energy of a soliton in 
d+1 dimensions is given by |C(~ \) [0], where £(s) is the zeta-function 

as) = — dtf-\K(t)-n ) (9) 







The effective action of an instanton on d dimensions is given by — 1£'(0). A proper time 
cutoff e on the lower limit of the integral can be introduced to regularise the theory |L7 . 
Inserting equation (H) gives 
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1 f°° — — 

C(-|) = / dkJ2(^(k) -5 7 (0)) + poles (10) 

7T J »y 
2 Z* 00 _ _ 

(>(0) = - dk k~ l y2(5Jk) - 5 JO)) + poles (11) 

7T JO „ 



The finite part is now in a form that can be evaluated numerically and the pole terms, which 
depend on the heat kernel coefficients, can be absorbed by renormalisation. 

Farhi et al. [[| have obtained expressions for the vacuum energy in 3 + 1 dimensions by 
subtracting the full first and second order phase shifts. This corresponds to subtracting an 
expansion of the heat kernel with nonlocal factors. Such an expansion does indeed exist [1 
although the local version is sufficient for removing the divergencies. 



III. CP(1) SOLITON 

The CP(1) model in 2 + 1 dimensions x^, fi = ... 2, has a complex field u with 
Lagrangian density 



1 d^ud^u 
g 2 (1 + mm) 2 



(12) 



where g~ 2 is an energy scale. The solutions to the field equations fall into families charac- 
terised by a topological index 0. The single soliton solution centered at the origin can be 
written u\ — a/z, where a is a constant and z is the complex coordinate x 1 + ix 2 . 
Perturbations £ about the soliton satisfy the equation 

-4D t D£ = k 2 Z (13) 

where D s = dg and D z = d z — 2u\{d z u{) / (1 + UiU\). There is one normalisable zero mode 
£ = z~ 2 and the remaining spectrum is positive. The first heat kernel coefficient for this 
operator can be evaluated from standard formulae [0,0 or by the Atiyah-Singer index 
theorem [[L3|] to be B\ — 1. 

The perturbation equation separates, and with the decomposition £ = (1 + a 2 / r 2 )u rn e im(l ', 

1 , (m + 2) 2 , , 2 



where the potential 



U m + ~ U 'm o U m + k U m = U (r)u m (14) 



4(m+l) 8a 2 

U H = 2 ^ 2 - / 2 , 2^2 15 

r z + a z yr z + err 



The boundary conditions are 



u _^ { ±J m+2 \ r -> (16) 

\ A(J m (kr) - Y m (kr) tan 5 m ) r -> oo '"' 



The potential does not meet the integrabihty requirements of Levinson's theorem |]TT[ and the 
numerical solution shows the unusual behaviour that S m (0) = tt for m > — 1 and 5 m (0) = —it 
for m < —2. 
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The first Born approximation to the phase shifts can be evaluted in terms of Bessel 
functions, 6$(k) = n(2m— ad a )I m {kot)K m {kot). The sum over m is not formally convergent, 
however the limit 



M 

lim E 8<£>{k)=* (17) 



m=—M 



is well defined and when substituted into equation (||) gives the correct conclusion that 
Bi = n Q . This definition of the m — > oo limit should also be used for actual phase shifts. 

The vacuum energy E can be obtained from equation (fLOl), with the observation that 
the first Born approximation cancels due to ([IT 



1 /"OO 00 

- dk E ( 6 rn(k) - 5 m (0)) (18) 

i7r J m =-oo 



where 2ir/g 2 is the classical contribution. The integrand is plotted in figure 1. After per- 
forming the integral numerically 

E = ^r- 0.248a" 1 (19) 

g 2 

The energy decreases as the width of the soliton decreases. However, if the width of the 
soliton changed with time, there would be an adiabatic solution of the form u = a(t)/z, but 
this solution has infinite action. The value of a must therefore remain fixed assymptotically, 
although it may decrease locally. 

The loop expansion parameter for the single soliton is effectively g 2 /ce, and so continuing 
to analysis to more than one loop would introduce further inverse powers of a [[TJJ. This 
raises the possibility that the energy has a minimum at some value of a. 
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FIG. 1. The integrand for the vacuum energy of the CP(1) soliton plotted as a function of 
k = ka, with angular modes up to M = 36. (The cusp at re = comes from the m = — 1 mode 
and, since <5'_i(0) ^ 0, 5_i(0) has not been subtracted.) 
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